Anisotropic Instabilities in Trapped Spinor Bose-Einstein Condensates 
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We theoretically investigate the effect of an anisotropic trap on the instability of the polar {mp = 
0) phase of a spin-1 Bose-Einstein condensate. By considering rigorously the spatial quantization, 
we show that the growth of the nascent ferromagnetic phase at short times becomes anisotropic 
with stronger oscillations in the magnetization correlation function along the unconfined direction. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) with active spin de- 
grees of freedom have proven to be a fertile ground for 
studying cooperative quantum many-body phenomena 
in condensed matter systems. Following the initial re- 
ahzation of BECs with spin-1 alkali atoms in an opti- 
cal trap [1] and description of the dynamics of the sys- 
tem [2, 3], theoretical and experimental investigations 
of spinor BECs have identified such diverse behavior as 
quantum phase transitions [4, 5], domain and topological 
defect formation [5-8] , coherent spin- mixing and amplifi- 
cation [9-11], and, most recently, dipolar effects [12-15]. 

In the experiment described in [5], a condensate of 
®^Rb atoms restricted to the ferromagnetic F = 1 hy- 
perfine manifold was prepared in the mp = state via 
the quadratic Zeeman effect. When the magnetic field 
was rapidly quenched, quantum fluctuations triggered 
the growth of unstable modes, leading to the formation 
of domains of magnetization with random orientations 
in the plane perpendicular to the spin-quantization axis. 
The instability is characterized by the inverse time con- 
stant F and corresponding length scale ^ = ^h/2mT, 
which are set by the strength of the spin-spin interaction 
[cf. Eq. (10)]. ^ coincides with the spin- healing length, 
with ^ = 2.4 /xm for the conditions of Ref. [5]. 

The exponential growth of the amplitude of a mode 
is proportional to e^^* and is controlled by the dimen- 
sionless gain parameter VL. In a homogeneous system, 
an unstable mode may be characterized by a wavevec- 
tor k. Neglecting the dipolar interaction, one finds [4] 
the gain parameter to be fJ = ^k'^(2 — fc^), where fc is 
measured in units of ^~^. Although the nature of the 
instability at early times (Fi ^ 4) has been successfully 
analyzed in the case of an infinite condensate [4, 6, 8], 
the effect of the trap geometry on the instability and 
spatial distribution of the associated magnetization has 
not been adequately addressed. This is a relevant ques- 
tion, since the experiment observed significant anisotropy 
in the magnetization correlation function which was at- 
tributed to the spatial anisotropy of the condensate. It 
was argued that as a result of the elongated shape of the 
condensate a small set of discrete modes in the trapped 
direction, rather than a continuum of states, led to pref- 
erential modulation of the magnetization in the confined 
direction. 



In this paper, we account rigorously for the spatial 
quantization effects in the development of the instability 
of the polar phase. Concentrating on the case of isotropic 
interactions within the F = \ manifold, we find that the 
correlation function [Eq. (15)] displays behavior precisely 
opposite to the one observed in experiment, i.e., displays 
a more prominent magnetization modulation in the un- 
confined direction. The inconsistency may be due to dis- 
tortions introduced by optical aberrations in the imaging 
device, and new experiments are underway [16]. 

To be specific, we consider a condensate whose thick- 
ness in the y-direction is comparable to or smaller than 
the spin-healing length so that the spin dynamics is ef- 
fectively constrained to the x — z plane. Moreover, we 
assume that the condensate is infinite in the z-direction 
but short in the x-direction due to a confining poten- 
tial. Note that the effect of a trap on single-component 
condensates has already been analyzed extensively [17]. 
However, the analysis of a condensate with a spinor or- 
der parameter is complicated by the fact that, in the case 
of isotropic interactions within the F = \ manifold, the 
SU{2) symmetry in spin space must be respected. In 
particular, the simplest approach in which we keep the 
functional form of r2(k) for an infinite system but restrict 
the allowed wavevectors to certain values by applying ei- 
ther hard- wall boundary conditions or periodic boundary 
conditions (PBCs) yields qualitatively incorrect results. 
In the case of hard-wall boundary conditions, kx takes 
on discrete values which exclude kx = 0, which results 
in il 7^ for all allowed values of k and contradicts the 
existence of a zero mode in the S'C/(2)-symmetric case 
[4]. Meanwhile, although the presence of the zero mode 
is preserved by choosing PBCs in the x-direction, this 
results in the wrong systematics of unstable modes (only 
an odd number of such modes arc possible, and the max- 
imal gain is identical for all modes with kx < 1). The er- 
rors introduced by either of these two simplifications are 
insignificant for wide traps {Lx 2> £,), but grow in impor- 
tance with decreasing Lx- For the typical trap widths 
of the order of a few 2tt^, these errors are substantial 
enough to affect the qualitative features of the instabil- 
ity development. As an illustration, we plot in Fig. 1 the 
gain parameter Q{kz) for a system with only a few modes 
spatially quantized in the x-direction by (a) PBCs and 
(b) a rigorous consideration of a harmonic trap. In the 
former case there are 2 (3, accounting for ±kx degener- 



acy) modes with equal maximum gain and another mode 
with a nearly-degenerate maximum. The exact degen- 
eracy of the different maxima for fc^, < 1 arises because 
the differences in k^ are compensated by different val- 
ues of kz satisfying k'^ = k^ + k^ = 1. In contrast, the 
latter is characterized by substantial differences between 
the maxima of Q for different modes. In time, this leads 
to an increasingly pronounced anisotropy of the magne- 
tization correlation function because a single mode with 
the largest gain dominates the growth of the nascent fer- 
romagnetic phase. Since this mode has a well-defined 
wavelength in the z-direction, increasingly strong oscil- 
lations between regions of positive and negative correla- 
tions are expected in the z-direction. Note that when 
Lx ^ Stt^ only a single mode corresponding to fc^ = is 
unstable and applying PBCs makes qualitatively correct 
predictions for the behavior of the magnetization corre- 
lation function. Indeed, this justifies our later treatment 
of the y-direction in a single-mode approximation, since 
Ly < 2tt^. 

In the following section we solve the equations of mo- 
tion for the amplitude of the confined mp = ±1 conden- 
sate during the exponential growth stage. On the basis 
of these solutions, in Sec. Ill we explore the anisotropy 
of the resulting magnetization correlation function. In 
doing so we emphasize the importance of the zero mode 
guaranteed by SU{2) symmetry, both to the spectrum of 
instabilities and to the stability of a spinor BEG in the 
absence of any spin-spin interaction. 



II. SHORT- TIME DYNAMICS 
A. Equations of motion 

In the absence of magnetic fields, an interacting spin-1 
Bose gas is described by the Hamiltonian [2] 
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where ^j? is a spinor whose components ^t^ create parti- 
cles in magnetic sublevel niF = +1, 0, or — I, h — ^^'^, 
F = iptF^ with F being the spin-1 matrices quantized 
along the z-axis, V{r) is the trap potential, and /u is 
the chemical potential. The interaction constants are 
related to the s-wave scattering lengths ap by cq = 
47rfi,^(ao -I- 202)73771, C2 = 47rfi,^(a2 — ao)/3m, and we 
assume that the spin-spin interaction is ferromagnetic 
(c2 < 0) as in ®'^Rb. To find the evolution of the ob- 
servables of the initially polar (tof = 0) state, we may 
impose two simplifications based on a zero-depletion ap- 
proximation whose domain of validity we consider later 
in the paper: 1) keep the mp = ±1 components to lowest 
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FIG. 1: (color online). Comparison of the gain parameter 
Q.{kz) for a two-dimensional spinor BEC in the x ~ z plane, 
for different models of confinement in the a;-direction. (a) 
Periodic boundary conditions (the simplest model) and (b) 
realistic confinement in a harmonic trap potential; Z/^/^ = 11 
in both cases. In (a) the solid and dashed curves correspond 
to fca; < 1 and kx > 1, respectively, and the two modes with 
n 7^ at fcz = are each doubly-degenerate. Note that all 
three modes with kx < 1 have the same maximum value of 
gain Q.. In (b), curves of the same parity (even for 1 and 3, 
odd for 2 and 4) show anti-crossing behavior due to mode- 
mixing, especially evident in the latter pair. The position 
of the arrow indicates the value of k^, corresponding to the 
maximum gain; the instability with this gain may become 
dominant before depletion becomes significant, see Sec. HID. 
The dashed curve in (b) is the fc^; = curve from (a), for 
comparison. 



order, i.e., to quadratic order in the Hamiltonian, and 2) 
treat ^0 as a classical real field '!/'o(r) satisfying the ordi- 
nary Gross-Pitaevskii equation for a scalar condensate: 
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With these assumptions the Heisenberg equations of mo- 
tion yield a closed system for ^+1, "^ _i: 

ihdt^+i = {Ho + C2n)*+i + C2n*li, (4) 

-ifiSt^Li = (Ho + C2n)'l'^ + C2?i^+i. (5) 
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Note that the latter form for Hq guarantees that ipo is an 
eigenstate oi Hq with zero eigenvalue. Furthermore, since 
"00 satisfying Eq. (3) is nodeless, the ground state energy 
and aU other eigenvalues of Hq are non-negative. As we 
will see later, this property is essential to respecting spin- 
conservation in this system in the absence of magnetic 
fields. 

Another important consequence is that a spinor BEC 
in the absence of spin-spin interactions is stable for any 
homogeneous spin state, an expected yet non-trivial fact. 
Although this is already clear from Eqs. (4), (5) with 
C2 = 0, it is worth seeing this explicitly from the multi- 
component Gross-Pitaevskii equation corresponding to 
Eq. (2) with ca = 0: 
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Because of the (spin-) rotational symmetry, the solution 
will be of the form ^o(r) = V'o(r)C7 where ^ is a position- 
independent spinor. This allows us to decompose small 
fluctuations about the mean-field as S'^ — (5^|| + S'^±, 
where S^h is along ( and *J(5*_l = (5*^*o = 0, so that 
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Here Hq is the same Hamiltonian from Eq. (7). The fact 
that Hq is bounded from below by zero therefore guaran- 
tees that in the absence of spin interaction, i.e., C2 = 0, 
the ground state energy of the BEC for any position- 
independent spin state (characterized by spinor C) is the 
same, and coincides with the one for a scalar BEC. 

Returning to the linearized quantum equations of mo- 
tion with C2 < 0, the density at the center of the conden- 
sate, no J defines a characteristic inverse time constant for 
the instability and corresponding length 



r = 
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for the spin-spin interaction. It will be convenient to 
measure times and distances in terms of T~^ and ^, re- 
spectively; we will use the notation 
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for the dimensionless quantities. We then arrive at the 
system of equations 



idA+i = {Hq - /?)$+! - p$li, 
-idr^li = {ho - p)<ILi - p<l+i 
where Hq — Hq/HT: 
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The observable quantity of interest is the transverse mag- 
netization correlation function, which we define using 

F± = Fx ± iFy as 



G^{r,r\t)^{F+{r,t)F-{r',t)) 
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where the dimensionless correlation function is given by 

9AQ,Q\r) = 2^p{Q)p{Q') (17) 

+ <l_i(e, r)<l+i(e', r) + <l_i(e, r)<l^i(e', r)) 

and expectation values are taken with respect to the ini- 
tial mp = ±1 vacuum state. In the following sections 
we derive the complete solution for g± from which the 
homogeneous solution for V = follows trivially, then 
specialize to the case where the condensate has a finite 
extent in one direction. 



B. General solution for the magnetization 
correlation function 



The solution of the coupled linear system in 
Eqs. (12), (13) can be written as 



$+i(e,T) = / dg' 
*Li(e,T)= [dg' 



Uig,g',T)^+i{g',0) + V*ig,g',T)'^l,ig',0) 

V{g,g\T)^+i{g\0) + U*{g,g\T)^l,{g\0) 
I 



(18) 
(19) 



with the initial conditions U{g, g' ,0) = 5(g — g'), thonormal basis (^a as 

V(q,q',0) = 0. We may expand U and V in an or- ^^, , , v-^ ^ \ / \ * ^ /\ /r>«\ 

^ ^ U{g,g',T) ^^^„{g)u,,p{T)^l{g'), (20) 

a, 13 

V{g, g', r)^Y. Vc.{qWp{t)^;{q') (21) 

a, 13 



where Wa/3(0) ~ Sap, 'Uq/3(0) — 0, while 

<l+i(e',o) = ^^^(e')a+i,7: (22) 
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$^(e',o) = ^^;(e')«^,,. (23) 



Hence we are primarily interested in computing Waj- 

Now, if we naively try to solve Eqs. (27), (28) by taking 
the time derivative of Eq. (28) and substituting Eq. (27), 
we get 



Then Eqs. (18), (19) become 
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^+i{q,t) = Y^ 
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</3a(e)«a/3(T)a+l,/3 + (^^ (£>X/3 (T)a^i ,3 

(25) 



Let Lpa be the cigcnstatcs of ho in Eq. (14) with eigenval- 
ues Ca- Substituting this expansion into Eqs. (12), (13) 
and introducing the linear combinations 



^ay — ^^a7 i '^ay, ^ay — ^07 ^07 

then gives 
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(26) 
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Here pafj ~ J dg fa{Q)p{Q)iPfi{Q) are the matrix ele- 
ments oi p(g). The corresponding initial conditions are 



Wa-yiO) = Za-yiO) = 6af 



(29) 



where Maf3 = £a^ai3 — '2'£aPa0 IS non-hcrmitian. Thus it is 
desirable to apply a transformation that renders the sys- 
tem hermitian. But we may note now that one mode of 
the system will be neither oscillating (positive eigenval- 
ues of M) nor growing (negative eigenvalues of M) . This 
follows from the earlier-noted existence of a mode wc la- 
bel as a = such that eo = 0, which implies M0/3 = 0. 
Since the first row contains only zeros, the determinant 
of M itself must be zero, which shows that zero is an 
eigenvalue of M. This is an important consequence of 
the SU{2) symmetry of the system [4]; specifically, the 
infinite- wavelength excitation creating an mp = ±1 pair 
is equivalent to a global rotation of the system, and hence 
does not cost any energy. 

By performing the similarity transformation w^-y = 

fa Waj, Zaj = £« Za-y for a ^ 0, WC find that the 
solution of Eqs. (27), (28) is governed by the hermitian 
matrix 



Mafj = elSafi - 2e]l^papey\ a, /3 ^ 0. 



(32) 



This is a particularly convenient formulation of the prob- 
lem, since the dimensionless magnetization correlation 
function in Eq. (17) is easily shown, with the help of 

Eqs. (24), (25) and the relations {amaO^jn'ii) = 5mm'5ap, 
to be 



9Aq,q',t) = 2^p{q)p{q') 
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(30) 



Let S be the unitary matrix that diagonalizes M, 
[5''''M5']aa' = E\5\x. Some of the eigenvalues E\ are 
negative, representing the unstable modes of the system. 
Since the oscillatory modes for which E\ is positive are 
quickly washed out by the exponential growth of the un- 
stable modes, we may restrict all summations over A to 
negative E\. For these we define Q,x = y/—E\, which 
we call the "gain" in reference to the connection with 
parametric amplifiers in quantum optics [18]. Then 
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In summary, the transverse magnetization correlation states ipa of Eq. (14) together with Eqs. (33)-(35). 
function is determined by Eq. (30) through the eigen- 



Note that previous results for the homogeneous case 
follow trivially from this formalism. li V = then n = 
ng is the solution to Eq. (3) and the eigenvalues and 
eigenstates of /ig [Eq. (14)] are ek = fc^ and (pk{g) oc e^^'^, 
respectively. Since p and M are diagonal, S is simply the 
identity matrix and the gain function is given by 



f^k= v/^2(2_fc2). 



(36) 



Eqs. (33)-(35) and Eq. (30) then give g±(Q,Q',T) = 
g±{Q — q',t) where 



ie^r) ex Y^ 
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(37) 



which is consistent with [6]. 

The formalism developed above is applicable to any 
trap geometry. In the case of tight confinement in one 
of the directions (the y-direction, for concreteness, where 
Ly ^ 27r^), one may use a single- mode approximation 
for the y-dependence of eigenf unctions tpa, as was done 
in previous works [8]. In the following sections we take 
into account the variation of (/?„ with x and z only, so 
that we may constrain the problem to the solution of 
differential equations in the x — z plane. 



III. ANISOTROPIC CONFINEMENT 
A. Motivation 

To gain an idea of what differences we should expect in 
the anisotropic case, note that for sufficiently large times 
we can approximate Eq. (37) as 
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Carrying out the sum using steepest-descent integration 
around the maximum gain at fc = 1 gives (in strictly two 
dimensions) 






(39) 



where Jq is a Bessel function. The dominant contribution 
to the correlation function is an average over a circle of 
wavevectors with unit magnitude, and thus manifests as 
radially decaying "oscillations" . 

As a qualitative first account of finite-size effects, 
we may consider quantizing the wavevectors in the x- 
direction as kx = 2Tm/lx while keeping kz continuous to 
simulate a highly anisotropic condensate; see Fig. 1(a). 
We can still carry out steepest-descent integration of 
Eq. (38) in the kz direction to obtain 
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cos(-\/l — k^z) 



ik^x 



(40) 



We recover Eq. (39) in the limit Z^: ^ oo by converting 
the sum over k^ to an integral. In contrast, for l^ ~ 27r 
only a few kx are allowed, and the correlation function 
develops some anisotropy. As noted in the introduction, 
however, for any /c^; < 1 the gain in Eq. (36) reaches the 
maximum il = 1 at some kz satisfying k'^ + k^ ~ 1 , which 
implies that several wavevectors are equally dominant in 
their contribution to the correlation function. Moreover, 
it is clear from the denominator in Eq. (40) that the result 
is highly dependent on the allowed values of kx^ i.e., the 
exact value of Ix- As we show in the following sections, 
both of these predictions are incorrect. To adequately 
account for the effect of confinement, we now consider 
more rigorously the spin dynamics in the presence of a 
harmonic potential. 



B. Spectrum of /lo 

Working in dimcnsionful units for a moment, suppose 
the condensate is infinite in the z-direction but has a 
finite extent Lx = 2_R in the x-dircction. We model 
the confining potential as V = }^mLiP'x^ ^ with the fre- 
quency determined by V(x = R) = ji = cqUq. Then in 
the Thomas- Fermi limit of vanishing kinetic contribution 
(oosc = y/h/muj ^ R), the solution to Eq. (3) is given 

by 



V'o 
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1/2 



(41) 



for |a;| < R and zero otherwise. Once |a:;| > R, the ex- 
act solution vanishes exponentially on the length scale 
S = {at^j2RY/^ < R (Ref. [17]). Correspondingly, the 
effective potential in the Hamiltonian Hq of Eq. (7) grows 
rapidly for |a:;| > R. This allows us to replace, in the limit 
(5 — > 0, the exact boundary conditions il){x -^ ±oo) — 
for the eigenfunctions of Hq with ^{x — ±R) — and 
simultaneously use Eq. (41) for -00 in Eq. (7): 
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The divergence of the second term in the Hamiltonian at 
X = ±R is regularized by the boundary condition ip{x = 
±i?) = 0. We also note that Eq. (3) can be solved directly 
by numerical integration [19] to give a solution that is 
valid for nonzero (5, but we opt for this simpler solution 
in order not to obscure the physics. 

Returning to dimcnsionlcss units, Eq. (14) becomes 



ho 



{xl — x^Y 



(43) 



where xq = R/S, and the solutions obey (p{x — ixg) = 0. 
The value of xq controls the dynamics in the trapped 
case, which is highlighted by the fact that it may be 



written as a simple combination of the three energy scales 
in the system: 
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The eigenstates and corresponding eigenenergies of this 
Hamiltonian are 

if>{Q) cx e*''->„(a;), eu^n ^ k^ + £„ (45) 

with ipn satisfying 

^n{x) = (46) 



dx'^ (xq — x'^Y 



and /_" dx |(p„(x)p = 1. If we make the transformation 



5„(x) = -1/1 2 <^" t^^^ ^ — ) ' 

V ^0 V ^0/ 



(47) 



then the eigenvalue problem reduces to finding K 
XqE, 0(??) such that 



-— r + K sech ?7 
diT^ 



Hv) = 



(48) 



and 0(?7 -^ ±00) = constant. This equation has been 
solved numerically; the first 4 eigenstates (converted back 
to (fnix)) and eigenvalues Kn are shown in Fig. 2 and Ta- 
ble I, respectively. Solutions to the square-well problem 
(width 2xo) obtained by neglecting the second term of 
Eq. (46) are also shown. The energy of the most rele- 
vant low-lying states differ significantly, illustrating the 
importance of accounting for the trap in a manner that 
properly preserves the zero mode. 

We have already seen the resulting gain function in 
Fig. 1(b), and it is worth noting that the presence of off- 
diagonal elements in M, Eq. (32), leads to mode-mixing 
and anti-crossing behavior among curves of the same par- 
ity. Such crossings appear in the nai've PBC quantiza- 
tion. There are remnants of the ±fca;-degeneracy from 
the homogeneous spectrum, especially near kz — 0. But 
most critical is the degree to which a single mode near 
/c^ ~ 1 is emphasized in the resulting gain function, so 
that we expect a greater dominance of the corresponding 
wavelength in oscillations of the magnetization correla- 
tion function. 



C. Magnetization correlation function 

Based on the gain parameter for the unstable modes 
[Fig. 1(b)], we claimed in the introduction the appear- 
ance of strong oscillations with a single wavelength in 
the z-direction. This is most dramatically illustrated by 
examining cross-sections of the correlation function at 




FIG. 2: First four eigenstates of Eq. (46) (solid lines), along 
with eigenstates of the square-well problem (dashed) . 
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4.62 


14.4 


29.1 


K7 
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7.40 


19.7 


37.0 



TABLE I: First four eigenvalues for of Eq. (46) compared to 
the square- well (sw) problem obtained by ignoring the second 
term; Kn = x^Sn is independent of the trap width. The 
bottom of the square well is taken to be at — 7r^/4 so that the 
lowest energy is zero as required by SU{2) symmetry. 



X — Q [Fig. 3(b)], where it is clear that oscillations of 
the magnetization correlation function in the ^-direction 
become stronger in time. The fact that this indicates 
the dominance of a single mode (contrary to that pre- 
dicted by PBCs) can be seen by plotting the positions of 
the first zero of (7j_(0, z,t) as shown in the lower curve 
of Fig. 4, which asymptotically approaches the value 
7r/2fc|, k* = 0.91. The latter corresponds to the maxi- 
mum of Vt{kz) [belonging to curve 4 in Fig. 1(b)]. We also 
show the first zero-crossing in the x-direction for compar- 
ison, although there is no meaningful wavelength in this 
direction and indeed the "oscillations" die out over time 
[Fig. 3(a)]. This is easily explained in the homogeneous 
picture where the dominant mode occurs at fe^ + fc^ = 1 • 
Since fc^ ~ 1, we expect the largest contribution in the 
x-direction to be k^ — 0, which has no oscillations. 

As we show in the next section, the underlying equa- 
tions of motion are valid only until t ~ 4. Meanwhile 
the asymptotic behavior is not realized until t ~ 15, as 
can be seen in Fig. 4. This can also be predicted on the 
basis of the gain curves in Fig. 1(b) by noting that the 
maximum gain of the top curve is separated from the 
maximum gain of the curve below by approximately 0.2; 
thus r = 15 represents roughly 3 time constants. Still, we 
expect even at fairly short times for the anisotropy to be 
visible in the full correlation function as shown in Fig. 5, 
where we plot the correlation function (7_l(0, q, t)/ ^J p{q) 
in the x — z plane. Note that, since we are interested in 
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FIG. 3: Plot of (a) g^{x,Q,T)/ ./p{x) and (b) g^{Q,z,T) di- 
vided by (/x (0, 0, r) . Arrows indicate the progression of curves 
with increasing r. Note the different scales on the horizontal 
axes. 
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FIG. 5: (color online). 




The density-normalized magnetiza- 

4 and 



tion correlation function g±(0, p,T)/^/p{g) at (a) r 
(b) T — 15. Dark regions represent alternating negative (red 
online) and positive (blue online) correlations, with strongest 
positive correlation at the center; the color scale is arbitrary. 
Large correlations were truncated to white to maximize con- 
trast in the rest of the plot. 



D. Depletion of the polar phase 



3.5 



FIG. 4: (color online). The position of the first node, de- 
fined by g±{0, zq,t) — for the z-direction (squares) and 
g±{xo, 0, r) = for the s-direction (circles), for the curves of 
Fig. 3. The dashed line indicates the asymptotic value, 1.7, 
predicted for the z-direction from the point of maximum gain 
in Fig. 1(b). 



the spin correlation, we have factored out y^. Smaller 
values of xq accelerate the development of anisotropy so 
that it may be conclusively discerned in a real experi- 
ment before depletion and nonlinear effects begin to play 
a significant role. 



In this section we consider for how long the underly- 
ing Eqs. (12), (13), which were based on a zero-depletion 
approximation, are valid. Neglecting thermal effects, the 
number of mp — particles taken out of the condensate, 
or equivalently the number of mp — ±1 pairs created, is 
given by iVpairs(T) = J dg npairs(£', t) where 

"pairsCe,^) = ($;^l(£.,T)l> + i(e,r)+l>Li(e,T)$_l(£>,r)> 

(49) 

= 2 ^ ^aig)y^*fiiQ)VcjiT)vl^{T). (50) 

Rather than evaluate what turns out to be a cumber- 
some expression, we note from Eqs. (26) and (30) that 
'T-pairsCff, t) < g^iQ^Q^T)! p{q). The expression is well- 
approximated by evaluating the integrand at ^ = 0, so 
that the condensate fraction k = Nps^i,-s/N is given by 



k{t) = 



g±(0,0,r) 



M) 



(51) 



where d is the dimensionality of the system. In the case 
of a quasi-2D trap, d = 2 and riQ—Jdy n(0, y, 0). The 
condition for the governing equations, and hence the re- 
sulting magnetization correlation function, to be a good 
description of the system is therefore k ^ 1. Since the 
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FIG. 6: The depletion function k,. Each unit of time repre- 
sents approximately 16 ms. We have chosen Ix ~ 11, h ~ 140, 
close to experimental conditions [5]. 
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mentally [5] , which poses the question of what is causing 
the observed preference for the trapped direction. An in- 
vestigation of the possible effects of optical aberration in 
the imaging device used in the experiment is underway, 
and the new experiments may shed light on the origin of 
the presently evident inconsistency [16]. 
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denominator in Eq. (51) represents the number of parti- 
cles in a spin correlation volume, physically this condition 
requires the magnetization per particle in each individ- 
ual domain to be small. Fig. 6 plots the depletion as 
a function of time, showing that the zero-depletion ap- 
proximation is quite good until r* ~ 4 (i.e., t* ~ 64 ms) 
for the conditions of [5] , after which the nonlinear effects 
not taken into account in the equations of motion are 
expected to play a significant role in the dynamics. The 
difference in gain Afi ~ 1/r* that can be resolved in the 
exponential growth of the different modes before deple- 
tion renders the underlying equations invalid is indicated 
by the length of the arrow in Fig. 1 (b) . 



IV. CONCLUSIONS 

We have analyzed the spin dynamics of a spatially 
anisotropic spinor BEG starting from the polar {mp — 0) 
state, which is unstable in the absence of a high magnetic 
field and develops ferromagnetic domains. We find that 
anisotropy should develop in the magnetization correla- 
tion function in the form of strong oscillations in the un- 
confined direction. We have ignored the effects of dipole- 
dipole interactions, which are inherently anisotropic but 
are an order-of-magnitude weaker than the spin-spin in- 
teractions [12]. Although dipolar effects can be impor- 
tant under certain circumstances [13, 15], preliminary 
investigation of dipolar forces in the system described 
here suggests that they will not significantly alter the 
magnetization correlation anisotropy. Nevertheless, our 
findings are contrary to what has been observed experi- 
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